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1.  Introduction 


The  US  Army  Research  Laboratory  (ARL)  is  investigating  the  advantages  of  group 
swarming  behavior  over  a  comparable  group  of  solitary  agents.  One  advantage  is 
the  ability  for  swarm  agents  to  localize  relative  to  the  group  using  spatial 
relationships  between  many  agents  to  achieve  accurate  relative  attitude  and  position 
information.1  This  is  particularly  important  in  GPS-denied  environments  where 
there  are  limited  options  for  absolute  positioning.2  In  these  cases,  it  is  possible  for 
a  swarm  with  relative  positioning  to  use  one  agent  with  absolute  localization  to 
grant  absolute  localization  to  the  entire  swarm.  Applications  also  exist  for  relative 
positioning  alone,  such  as  collision  avoidance,3  formation  flying,4  and  patterned 
weapon  delivery.5 

One  method  to  achieve  relative  localization  is  radio  direction  finding  (RDF).  There 
are  several  ways  to  perform  RDF,  but  commonly  the  phase  differences  between  the 
elements  of  an  antenna  array  are  used  to  determine  a  radio  signal’s  angle  of  arrival 
(AoA).6  An  agent  can  then  determine  its  relative  attitude  and  position  using  known 
angles  to  the  other  swarm  agents.  RDF  has  been  used  since  World  War  I7  and  has 
many  applications  such  as  ship  and  aircraft  navigation,  search  and  rescue,  wildlife 
tracking,  location  of  illegal  radio  transmissions,  missile  guidance,  radar  systems, 
and  antiradiation  missiles.8 

Despite  RDF’s  long  history  and  many  successful  applications,  a  suitable 
commercial  off-the-shelf  device  for  swarms  of  small  agents  does  not  exist.  Larger 
systems  for  search  and  rescue9  and  radar  cannot  be  easily  adapted  to  small  swarm 
agents,  while  smaller  automotive  radar  systems  have  limited  range.10  Instead  of 
using  the  phase  interferometry  approach,  many  systems  use  time  difference  of 
arrival  or  received-signal-strength-indication  methods  that  do  not  achieve  high 
accuracy.1 1  Phase  interferometry  using  antenna  arrays  can  be  difficult  to  implement 
because  antenna  characteristics  and  multipath  effects  can  significantly  degrade 
performance.  Despite  these  problems,  literature  suggests  that  a  reliable  RDF  system 
could  be  designed  for  swarm-relative  localization.12  There  are  2  significant  reasons 
why  many  swarm  applications  present  fewer  design  challenges  than  other  RDF 
applications.  First,  high-altitude  swarms,  such  as  swarms  of  unmanned  aerial 
vehicles  (UAVs)  or  munitions,  will  not  be  significantly  affected  by  multipath. 
Secondly,  many  swarms  can  have  antenna  arrays  attached  to  the  agents  in  a 
repeatable  manner,  simplifying  antenna  characterization  issues.  Compare  this  with 
a  device  such  as  a  cell  phone,  whose  possible  orientation  and  proximity  to  other 
objects  makes  assumptions  about  antenna  performance  problematic. 
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This  report  shows  the  feasibility  of  RDF  for  swarm-relative  localization  by 
demonstrating  RDF  with  a  simple  2-element  antenna  array  with  a  small  standalone 
software-defined  radio  (SDR).  First,  the  theory  of  phase-interferometry  RDF  is 
presented.  Next,  an  RDF  system  is  described  and  simulated,  including  the  practical 
considerations  of  the  radio  receiver.  The  software  design  and  hardware  setup  of  an 
RDF  laboratory  experiment  is  explained  next,  followed  by  the  experimental  results. 
Finally,  the  results  are  summarized  and  future  research  directions  are  explored. 

2.  Phase-Interferometry  Radio  Direction  Finding 

The  AoA  of  an  RF  signal  can  be  determined  using  the  system  in  Fig.  1.  A1  and  A2 
indicate  the  position  of  2  antennas  spaced  a  distance  d  apart.  The  red  lines  indicate 
the  RF  signal’s  direction  of  propagation  from  the  source  to  the  antennas.  If  the 
distance  from  the  signal  source  to  the  antennas  is  much  greater  than  d,  these  lines 
can  be  assumed  to  be  parallel.  The  AoA,  9,  is  shown  as  the  angle  off-center  of  the 
incoming  RF  signal  that  intersects  A2.  Since  9  +  a  —  90°  and  9  +  aPA1A2  —  90°, 
then  aPA1A2  is  also  9.  The  length  of  PA2  is  the  difference  in  distance  that  the  RF 
signal  must  travel  to  A2  as  compared  with  A1.  Given  an  RF  signal  with  wavelength 
A,  assuming  d  <  A/2  and  —90°  <  9  <  90°,  the  phase  difference  between  the 
received  signals  at  A1  and  A2  is  given  by 


(1) 


where  —  n  <  A <p  <  n.  Here  and  in  the  rest  of  the  report,  phases  are  in  radians,  while 
angles  are  in  degrees.  Given  a  measured  phase  difference  between  A1  and  A2,  the 
angle  of  arrive  is  then 


(2) 


Fig.  1  AoA  system  with  2  antennas 


Approved  for  public  release;  distribution  is  unlimited. 


2 


From  Eq.  1,  if  d  >  A/2,  A <p  can  exceed  the  range  of  —n  <  A <fi  <  n.  This  causes  an 
ambiguity  problem  since  values  of  A <p  and  A <p  +  2n  cannot  be  distinguished  from 
each  other  at  the  antennas,  leading  to  multiple  solutions  for  9.  To  correct  this 
problem,  an  integer  /  can  be  added  to  Eq.  2  to  account  for  the  phase  roll-overs: 


9  —  sin  1 


(3) 


For  a  range  of  —9max  <  9  <  9max ,  the  maximum  distance  of  PA2  is  ±d  sin  9max. 
Examining  increasing  values  of  9 ,  the  first-phase  roll-over  will  occur  when  A <p  = 
7r,  corresponding  to  PA2  —  2/2.  Each  additional  roll-over  will  occur  at  an 
additional  phase  change  of  2n,  corresponding  to  additional  multiples  of  A  in  length 
PA2 .  Therefore,  for  a  range  of  -9max  <9  <  9max,  the  range  of  /  will  be  -Imax  < 
I  <  imax ,  with  Imax  given  as 

j  _  I  d  sin  9max+A./2  .. 

^max  —  ^  >  v^v 

with  [xj  indicating  the  floor  of  x.  A  typical  method  to  determine  /  is  to  use  3 
antennas,  with  distances  d12  and  d13  denoting  the  spacing  between  the  first  and 
second  and  first  and  third  antennas,  respectively,  leading  to  2  phase  measurements, 
A <p12  and  A013.13  Rearranging  Eq.  3  for  each  pair  of  antennas  gives 


A012-27r/12  =  27r(^f^) 


(5) 


and 


A0i3  2tt/13 


(6) 


Solving  for  A 013  in  terms  of  A <p12  gives 


A013=-^A012--^27t/12  +  27t/13. 

ai2  a  12 


(7) 


In  many  AoA  applications,  the  direction  finding  system  has  no  control  over  the 
signal  being  detected,  but  in  applications  where  one  has  control  over  the  transmitted 
signal,  an  alternative  way  of  resolving  the  ambiguity  presents  itself.  Instead  of  using 
2  different  antenna  spacings,  2  different  transmit  frequencies  can  be  used. 
Equations  5  and  6  now  become 

A01-27t/1  =  27t(^)  (8) 

and 

A02-2tt/2  =2tt(^),  (9) 

and  Eq.  7  becomes 

A02  =  ^A0X  +  2nl2.  (10) 

^2  'L 2 
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Equation  10  indicates  a  linear  relationship  between  the  phase  measurements  A <p1 
and  A02  f°r  each  possible  combination  of  7X  and  72,  which  will  be  referred  to  as 
phase  lines.  An  example  scenario  showing  values  of  A<p1  and  A<p2  for  —80°  <  9  < 
80°  is  shown  in  Fig.  2  for  d  —  26.6  cm,  Xx  —  12.5  cm,  and  X2  —  15.625  cm.  In 
this  case,  both  7X  and  /2  are  in  the  range  —2  <  /x,  /2  <  2  for  —80°  <  9  <  80°, 
resulting  in  the  phase  lines  shown  in  Fig.  3.  A  total  of  100  phase  measurements 
from  Fig.  2  are  plotted  as  (A<p1,A<p2)  phase  points  that  occur  on  their  respective 
phase  lines.  The  line  of  a  phase  point  determines  its  values  of  l1  and  72,  which  can 
then  be  inserted  into  Eqs.  8  and  9  to  unwrap  the  phase  measurements  and  determine 
an  unambiguous  value  for  9.  Figure  4  shows  this  process  by  displaying  the  original 
phase  measurements,  the  correction  factor  of  —2nl  from  Eqs.  8  and  9,  and  the 
unwrapped  phases  for  both  Xx  and  X2  found  by  adding  the  correction  factor  to  the 
phase  measurements. 


Fig.  2  Example  phase  measurements  for  2  wavelengths 


A 
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Fig.  3  Example  phase  line  plot  with  phase  points 
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Fig.  4  Unwrapping  phase  measurements  for  the  first  wavelength  (top)  and  second 
wavelength  (bottom) 
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This  example  showed  an  ideal  case  with  perfect  measurements.  Figures  5-8  show 
the  same  example  with  JV (0,0.025)  Gaussian  noise  added  to  the  A <p  measurements. 
In  Fig.  6  the  phase  points  (A 01;  A <p2)  no  longer  lie  directly  on  the  phase  lines.  This 
means  that  we  can  no  longer  directly  determine  the  values  of  It  and  I2  by  solving 
Eq.  10  (i.e.,  matching  the  phase  point  to  its  respective  phase  line).  Instead,  the  phase 
line  closest  to  the  phase  point  is  used  to  determine  its  values  of  lt  and  I2.  The 
distance  5  from  point  (A 01(  A <p2)  to  the  phase  lines  in  Eq.  10  is 


(11) 


Using  5  to  determine  the  most-likely  values  of  I{  and  I2  for  each  phase  point, 
Fig.  7  shows  the  correction  factor  of  —2nl  and  the  unwrapped  phases  for  A1  and  A2 
measurements  with  the  final  angle  calculations  and  angle  error  displayed  in  Fig.  8. 
The  angle  error  increases  with  9  because  for  larger  angles  a  small  change  in  phase 
corresponds  to  a  large  change  in  6 . 
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Fig.  5  Noisy  phase  measurements  for  2  wavelengths 
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Fig.  6  Phase  lines  with  noisy  phase  points 


20 


A  0, 

-2  1 

-  7T  i 

Unwrapped 

A</>i 

Q  (Degrees) 


20 


-  A()>2 

-2  I 

-  7 T  2 

Unwrapped 

A  (^)2 

Q  (Degrees) 

Fig.  7  Unwrapping  noisy  phase  measurements  for  the  first  wavelength  (top)  and  second 
wavelength  (bottom) 


Approved  for  public  release;  distribution  is  unlimited. 

7 


True  q  (Degrees) 


True  q  (Degrees) 

Fig.  8  Recovered  angle  vs.  the  true  angle  (top)  and  the  recovered  angle  error  (bottom)  for 
noisy  phase  measurements 

As  the  level  of  noise  increases,  phase  points  may  occur  closer  to  the  wrong  phase 
line,  resulting  in  erroneous  values  for  /x  and  72.  By  designing  the  system  to 
maximize  the  spacing  between  the  phase  lines,  the  probability  of  these  errors  can 
be  minimized.14  The  spacing  between  2  phase  lines,  p,  with  integer  phase 
corrections  7X1  and  /12  for  the  first  line  and  /21  and  /22  for  the  second  line,  is  given 
by 

I  jt27r(/2i-/n)  +  27r(/12-/22)| 

P  =  - - 1— - '  (12) 


The  upper  plot  of  Fig.  9  shows  the  minimum  phase  line  spacing  for  varying  values 
of  A2  while  A1  is  held  constant  at  12.5  cm.  The  bottom  plot  shows  the  maximum 
values  of  It  and  /2  used  to  calculate  the  minimum  phase  spacing  for  9max  —  80°. 
Some  of  the  discontinuities  in  the  upper  plot  are  caused  by  changes  in  the  maximum 
values  of  7X  and  72  shown  in  the  bottom  plot.  For  some  values  of  A2,  the  phase 
spacing  is  zero,  indicating  that  some  of  the  phase  lines  overlap,  making  it 
impossible  to  resolve  phase  ambiguities.  Previous  examples  used  A2  —  15.625  cm, 
which  is  indicated  in  Fig.  9  to  be  a  local  maximum,  creating  relatively  wide  phase¬ 
line  spacings.  The  case  of  A2  —  11.538  cm,  also  marked  in  Fig.  9,  is  an  example 
with  relatively  small  phase  spacing,  resulting  in  the  phase  lines  in  Fig.  10. 
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Fig.  9  Minimum  phase-line  spacing  vs.  wavelength  (top)  and  the  corresponding  maximum 
integer  ambiguity  values  (bottom) 
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Fig.  10  Example  of  nonideal  phase-line  spacing 
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Although  in  many  direction-finding  applications  the  assumption  that  the  distance 
to  the  radio  source  is  much  greater  than  d,  in  swarming  scenarios  it  is  possible  that 
agents  may  be  quite  close  to  each  other.  Therefore,  it  is  necessary  to  explore  the 
accuracy  of  the  previous  assumption  that  the  incoming  radio  signals  to  both 
antennas  are  parallel.  Figure  11  shows  the  geometry  of  an  RDF  system  where 
incoming  RF  signals  are  not  assumed  to  be  parallel.  As  in  Fig.  1,  Ax  and  A2  are  the 
positions  of  the  2  antennas.  6  is  the  AoA  but  is  now  measured  from  the  midpoint 
between  Ax  and  A2  labeled  M.  rl.  rm,  and  r2  are  the  distances  from  the  RF  source 
to  At,  M,  and  A2,  respectively.  Using  the  law  of  cosines, 

ri  =  Jr£+^-rmdcos/3  (13) 

and 


r2  = 


d2 

+  —  —  rmd  cos  a. 


where  /?  =  90°  —  9  and  a  —  90°  +  6. 


Fig.  11  AoA  system  with  source  close  to  antenna  array 


(14) 


The  top  plot  of  Fig.  12  shows  the  error  in  the  calculated  AoA  using  the  parallel 
assumption  from  Fig.  compared  with  the  more  accurate  geometry  in  Fig.  1 1  for  A  — 
12.5  cm,  d  —  A,  and  rm  —  10  m.  The  bottom  logarithmic  plot  shows  the  maximum 
error  versus  the  source  range  rm  for  3  different  choices  for  antenna  spacing  d.  At  a 
distance  of  0.5  m,  even  the  worst  case  has  an  error  of  less  than  1°,  showing  that 
even  for  relatively  close  distances  the  parallel  assumption  is  very  accurate. 
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Fig.  12  AoA  error  of  the  parallel  assumption  for  an  example  case  (top)  and  the  maximum 
error  vs.  rm  for  3  values  of  d  (bottom) 

3.  Experiment  Setup 

Now  that  RDF  theory  has  been  presented,  a  practical  experiment  to  test  the 
performance  of  a  RDF  system  is  described  in  the  following.  Figure  13  shows  the 
experiment  setup.  An  RF  signal  generator  is  used  to  produce  the  source  signal, 
which  is  transmitted  by  a  horn  antenna  in  an  anechoic  chamber.  The  receiving 
antenna  array  and  SDR  are  mounted  on  a  rotary  actuator.  Figure  14  shows  a  picture 
of  the  antenna  array  mounted  on  the  actuator  in  the  anechoic  chamber,  and  Fig.  15 
shows  a  close-up  of  the  antenna  array.  A  PC  controls  the  actuator  to  vary  the  AoA 
and  records  data  from  the  SDR  for  RDF  postprocessing.  The  SDR  is  discussed, 
followed  by  an  overview  of  the  control  software. 
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Aiiechoic  Chamber 


Fig.  13  RDF  experiment  setup 


Fig.  14  Receiver  setup  in  the  anechoic  chamber 
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Fig.  15  Close-up  of  an  example  antenna  array.  The  antenna  in  the  middle  is  used  for  a 
calibration  process  described  later. 

SDRs  are  versatile  platforms  that  are  ideal  for  custom  applications  such  as  RDF. 
Whereas  traditional  radio  systems  are  entirely  implemented  in  hardware,  SDRs 
digitize  the  RF  signals  and  then  process  them  in  software.  This  makes  it  possible  to 
directly  determine  the  phase  differences  necessary  for  AoA  determination.  Building 
on  ARL’s  previous  experience  with  Universal  Software  Radio  Peripheral  (USRP) 
SDRs,1516  the  USRP  E310  shown  in  Fig.  16  was  chosen  for  this  application.17  The 
E310  can  act  as  a  standalone  Linux  system  with  2  receivers  and  2  transmitters 
operating  between  70  MHz  and  6  GHz.  Its  small  size  and  standalone  operation 
make  it  a  potential  choice  for  a  fielded  UAV  RDF  system. 
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Fig.  16  USRP  E310 

Figure  17  shows  a  simplified  block  diagram  of  the  E310.17  Received  RF  signals 
travel  through  filter  banks  before  entering  the  Analog  Devices  9361  RF  integrated 
circuit  (RFIC).  Inside  the  RFIC,  the  signal  is  amplified,  Inphase/Quadrature  (IQ) 
demodulated  to  an  intermediate  frequency,  and  digitally  sampled.  The  samples  are 
passed  to  the  PL  (programmable  logic)  on  the  Zynq  integrated  circuit,  where  they 
are  DDC  (digitally  down  converted)  to  baseband  and  decimated.  The  final  IQ  data 
is  then  passed  to  the  Zynq  PS  (processing  system)  running  OpenEmbedded  Linux 
for  further  processing. 


Fig.  17  USRP  E310  block  diagram 


Although  the  IQ  demodulation  in  the  E310  involves  a  number  of  steps,  the  basic 
functionality  is  shown  in  Fig.  18. 
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Fig.  18  IQ  demodulator  block  diagram 


The  input  RF  signal  x(t)  with  carrier  frequency  fc  can  be  represented  using  IQ 
components  as 

x(t)  =  /(t)  cos(27r/ct)  —  Q(t)  sin(27r/ct).  (15) 

Starting  with  the  demodulation  of  the  /  data,  Eq.  15  gives 

=  2cos(2nfct)(I(t )  cos(2nfct )  —  Q(t)  sin(27r/ct)).  (16) 

Using  the  double  angle  formulas 

cos(20)+l 

cos2(0)  = - - -  (17) 

and 

sin(0)cos(0)  =  sin(20),  (18) 

hj  (t)  becomes 

/i/(t)  =  /(t)  +  I{t)cos{4nfct )  —  Q(t)  sin(47r/ct).  (19) 

This  leaves  the  baseband  I  data  with  additional  signals  at  twice  the  carrier 
frequency.  Filtering  out  these  higher  frequencies  using  a  low-pass  filter  will  restore 
the  original  /  data.  Similarly,  for  the  Q  data 

hQ(t)  —  —2sin(2nfct)(I(t)  cos(2nfct )  —  Q(t)  s\n(2nfct)').  (20) 

With  the  additional  double  angle  formula 

sin2(6>)  =  1~C02s(2e),  (21) 

hQ  ( t )  becomes 

hq{t)  —  —I(t)sin(4nfct)  +  Q(t)  —  Q(t )  cos(47r/ct).  (22) 

Here  the  baseband  Q  data  is  left  with  additional  high-frequency  signals.  Filtering 
out  these  high  frequencies  will  restore  the  original  Q  data. 
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IQ  demodulation  can  be  applied  to  the  signals  in  this  RDF  experiment.  A  single 
frequency,  fc  +  fb,  is  used  as  the  RF  source.  Through  IQ  demodulation,  it  is  down- 
converted  to  a  fb  baseband  frequency  while  the  phase  information  is  preserved.  To 
understand  this  process,  we  convert  the  source  signal  to  the  IQ  format  in  Eq.  15. 
Let  the  source  be  represented  as 

x(t)  =  cos(27rt(/c  +  fb )).  (23) 

The  received  signals  at  antennas  Ax  and  A2  are  then 

*i(  0  =  cos(27rt(/c  +  fb)  +  0i)  (24) 

and 

x2(t)  =  cos(2?rt(/c  +  fb)  +  02).  (25) 

Using  the  cosine  sum  formula 

cos(ct  +  /?)  =  cos  a  cos  (3  —  sin  a  sin  (3,  (26) 

with  a  —  2ntfb  +  0  and  /?  =  2ntfc,  Eqs.  24  and  25  can  be  written  as 

xx(t)  =  cos(2ntfb  +  0i  )  cos(2ntfc  )  —  sin(2ntfb  +  0X)  sin(27rt/c)  (27) 

and 

x2(t)  =  cos(2ntfb  +  02  )  cos(2ntfc  )  —  sin(27rt/fe  +  02)  sin(27rt/c).  (28) 

This  is  the  IQ  signal  representation  in  Eq.  15,  with  /(t)  =  cos(2ntfb  +  0)  and 
Q(t )  =  sin(2ntfb  +  0).  The  phase  is  then 

0  =  tan-1f|.  (29) 

Thus,  the  recovered  IQ  data  preserves  the  phase  information  of  the  antenna  signals 
and  can  be  used  to  calculate  the  AoA. 

To  acquire  the  IQ  data,  the  USRP  E310  was  programmed  using  GNU  Radio18  to 
configure  its  2  receivers  and  send  IQ  data  over  user  datagram  protocol  (UDP)  to 
the  control  PC.  The  program,  test_udp.py,  is  included  in  Appendix  A.  PuTTY,19  an 
open-source  terminal  emulator  on  the  control  PC,  was  used  to  interface  with  the 
E310  and  run  the  GNU  Radio  programs.  A  Lab  VIEW20  program,  shown  in  Fig.  19, 
was  developed  to  receive,  analyze,  display,  and  save  the  IQ  data  over  UDP  from 
the  E310.  It  also  controlled  the  rotary  actuator  by  executing  commands  through  an 
emulated  serial  port.  AoA  experiments  were  automated  by  configuring  the  actuator 
to  sequence  through  a  series  of  angles  while  saving  data  from  the  E310.  Although 
the  Lab  VIEW  program  performed  some  basic  processing  during  the  experiments, 
the  final  analysis  was  performed  on  the  recorded  data  using  Matlab.21 
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Fig.  19  Example  screenshot  of  the  Lab  VIEW  program  used  to  control  the  AoA  experiment 

4.  Experiment  Results 

Several  experiments  were  performed  to  characterize  the  performance  of  this  RDF 
system.  The  first  was  a  high-resolution  experiment  using  standard  antenna  spacing 
and  power  levels.  The  second  tracked  performance  as  the  power  level  was 
decreased.  The  third  experiment  investigated  the  effect  of  varying  the  antenna 
spacing.  The  last  experiment  explored  the  feasibility  of  determining  the  SDR’s 
unknown  initial  phase  offset. 

In  the  first  experiment,  2  wavelengths  were  used  to  resolve  phase  ambiguities  using 
the  parameters  listed  in  Table  1.  Figure  20  shows  example  IQ  data  from  the  E310 
for  the  2  antennas  at  a  single  8.  Phases  were  extracted  from  the  IQ  data  using 
Eq.  29  and  shown  in  the  top  plot  of  Fig.  21.  The  middle  plot  shows  the  difference 
between  the  phases.  Values  outside  of  ±n  occur  due  to  phase  rollovers.  The  bottom 
plot  shows  the  corrected  phase  difference  with  all  of  the  values  wrapped  inside  of 
+7T.  The  average  of  these  phase  differences  is  used  for  AoA  calculations.  At  a 
modest  sampling  rate  of  50  kHz,  27,433  samples  were  averaged  over  a  period  of 
0.55  s  to  produce  the  final  A <p  measurement  for  a  given  8.  Figure  22  shows  the 
calculated  A <p  measurements  for  —45°  <  8  <  45°  with  steps  of  2.5°,  making  a  total 
of  36  A <p  samples  for  each  wavelength.  From  these  measurements,  the  phase  points 
were  plotted  with  their  phase  lines  in  Fig.  23.  As  in  Fig.  6,  these  points  are  not 
positioned  exactly  on  their  phase  lines,  so  the  closest  phase  line  was  determined 
using  Eq.  1 1  to  resolve  the  integer  ambiguities  I±  and  I2 .  These  values  are  then  used 
to  unwrap  A<p1  and  A<p2  as  shown  in  Fig.  24,  and  converted  to  the  estimates  of  8  in 
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Fig.  25.  The  results  were  fairly  accurate,  with  a  maximum  error  of  about  5°.  More 
encouraging,  however,  are  the  results  in  Fig.  26  showing  the  repeatability  of  the 
same  experiment  performed  3  times,  which  had  an  average  variance  of  0.008°  over 
the  3  trials  of  the  36  angle  estimates  in  both  X1  and  X2 .  This  indicates  that  significant 
increases  in  RDF  accuracy  may  be  achievable  through  calibration.  It  is  difficult  to 
determine  the  effectiveness  of  a  calibrations  procedure,  however,  until  additional 
experiments  are  performed  in  more-realistic  outdoor  environments. 


Table  1  Experiment  1  parameters 


Parameter 

Value 

Definition 

12.5  cm 

Wavelength  1 

^2 

15.625  cm 

W avelength  2 

d 

26.6  cm 

Antenna  spacing 

Grx 

-3  dBm 

Signal  generator  gain 

Grxi 

70  dBm 

Receiver  1  gain 

Grx  2 

70  dBm 

Receiver  2  gain 

fs 

50  kHz 

Sampling  frequency 

h 

1  kHz 

Demodulated  baseband  frequency 

0.4 


Time  (s) 


0.4 


Time  (s) 


Fig.  20  IQ  signals  from  antenna  1  (top)  and  antenna  2  (bottom) 
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Fig.  21  Phase  calculated  from  the  IQ  data  (top),  phase  difference  (middle),  and  corrected 
phase  difference  (bottom) 


Q  (Degrees) 


Fig.  22  Experimental  phase  measurements  for  2  wavelengths 
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Fig.  23  Phase  lines  with  experimental  phase  points 
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Fig.  24  Unwrapping  experimental  phase  data  from  first  wavelength  (top)  and  second 
wavelength  (bottom) 
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Fig.  25  AoA  experiment  results  using  2  wavelengths  (top)  and  AoA  error  (bottom) 
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Fig.  26  Repeatability  of  experiment  results  using  XI  (top)  and  X2  (bottom) 
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The  experiment  used  the  phase  lines  closest  to  the  phase  points  to  determine  the 
values  of  7X  and  72,  as  shown  in  Fig.  23.  In  the  following  experiment,  phase  error 
caused  some  of  the  phase  points  to  appear  closer  to  the  wrong  phase  lines,  resulting 
in  erroneous  values  of  I±  and  /2.  Here,  12.5  cm  was  used  for  A1  and  14.58  cm  for 
A2.  Figure  27  shows  the  phase  measurements  and  Fig.  28  shows  the  corresponding 
phase  points,  with  the  measurements  that  will  result  in  erroneous  integer  ambiguity 
resolution  marked.  In  Fig.  28  these  erroneous  points  are  closest  to  the  phase  line 
corresponding  to  7X  =  —  2  and  72  =  —1.  This  leads  to  the  unwrapped  phase 
differences  in  Fig.  29,  with  the  erroneous  —2nl  correction  factors  indicated.  This 
creates  obvious  phase  difference  errors  in  the  unwrapped  phases  and  angle  errors 
in  Fig.  30. 


A0, 
A  (f>2 


Q  (Degrees) 


Fig.  27  Experiment  phase  measurements  with  erroneous  points  marked 
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Fig.  28  Phase  lines  and  phase  points  with  erroneous  points  marked 
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Fig.  29  Experimental  results  using  2  wavelengths  with  erroneous  points  marked 
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Fig.  30  AoA  experiment  results  using  2  wavelengths  (top)  and  AoA  error  (bottom)  with 
erroneous  points  included 


These  types  of  errors  can  be  corrected  using  a  threshold  for  the  maximum  difference 
between  one  A <p  measurement  and  the  next.  In  this  case,  the  erroneous  integer 
ambiguity  resolution  caused  a  jump  in  the  unwrapped  phase  difference  of  about  2n. 
Using  a  threshold  of  A<pmax  —  n,  if  the  difference  between  the  previous  A <p 
measurement  and  current  A <p  measurement  is  greater  than  A <pmax,  the  closest  phase 
line  is  discarded,  and  the  next  closest  phase  line  is  used  to  resolve  the  integer 
ambiguity.  This  additional  step  results  in  the  corrected  unwrapped  A <p 
measurements  in  Fig.  31  and  9  measurements  in  Fig.  32.  In  a  field  application,  the 
higher  the  AoA  sampling  rate,  the  lower  the  difference  between  successive  A <p 
measurements.  Thus,  using  high  A <p  sampling  rates  could  make  this  type  of  error 
correction  very  robust  due  to  the  tight  bounds  on  A<pmax.  It  may  also  be  possible  to 
use  inertial  measurements  from  the  transmitter  and  receiver  agents  to  estimate  the 
difference  in  A <p  between  measurements,  providing  additional  bounds  on  A 4>max. 
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Fig.  31  Corrected  phase  unwrapping  of  first  wavelength  (top)  and  second  wavelength 
(bottom) 


Fig.  32  Corrected  AoA  results  (top)  and  AoA  error  (bottom)  for  2  wavelengths 
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The  next  experiment  determined  the  effects  of  transmit  power  on  AoA  error. 
Figures  33  and  Fig.  34,  corresponding  to  Figs.  20  and  Fig.  21,  show  example  noisy 
IQ  and  phase  data.  The  spectrum  of  this  noisy  data  is  shown  in  the  top  plot  of 
Fig.  35.  Here,  fb  is  around  2  kHz  and  can  be  clearly  seen  as  the  frequency  with  the 
highest  amplitude  in  the  plot.  Since  the  signal  is  so  narrowband,  a  narrow-bandpass 
filter  can  be  used  to  filter  out  most  of  the  noise  in  the  rest  of  the  spectrum.  The 
bottom  plot  of  Fig.  35  shows  the  frequency  response  of  a  16-pole  Butterworth  filter 
with  a  200-Hz  passband  designed  to  filter  the  IQ  data.  The  filtered  IQ  data  are 
shown  in  Fig.  36  with  corresponding  phase  signals  in  Fig.  37.  Clearly,  filtering  has 
greatly  improved  the  phase  measurements.  Using  this  filtering  method,  AoA 
experiments  were  performed  at  various  transmit  power  levels.  The  top  plot  of 
Fig.  38  shows  the  AoA  error  of  4  experiments,  2  using  a  transmit  power  of 
-53  dBm  and  2  using  -73  dBm.  Even  at  -53  dBm,  the  error  is  relatively  low  and 
appears  to  be  repeatable.  At  -73  dBm,  the  error  is  significantly  higher  but  still 
somewhat  repeatable.  It  would  be  expected  that  as  the  noise  increases,  error  caused 
by  noise  would  overtake  other  error  sources,  resulting  in  random,  nonrepeatable 
AoA  errors.  The  plot,  however,  unexpectedly  shows  somewhat  repeatable  error 
even  in  the  -73-dBm  case.  This  interesting  phenomenon  requires  further  study.  The 
bottom  plot  of  Fig.  38  shows  the  average  error  over  a  range  of  transmit  powers. 
Using  information  about  the  test  setup,  it  is  possible  to  convert  this  transmit  power 
to  an  estimated  distance. 
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Fig.  33  IQ  data  from  antenna  1  (top)  and  antenna  2  (bottom)  for  a  transmit  power  of 
-53  dBm 
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Fig.  34  Phase  calculated  from  noisy  IQ  data  (top),  phase  difference  (middle),  and  corrected 
phase  difference  (bottom) 


-50 


Q.  -150  _ | _ | _ | _ | _ 

<  0  5  10  15  20  25 

Frequency  (kHz) 


Frequency  (kHz) 

Fig.  35  Spectrum  of  example  noisy  IQ  data  (top)  and  frequency  response  of  a  filter  designed 
to  clean  the  IQ  data  (bottom) 
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Fig.  36  Filtered  IQ  data  from  antenna  1  (top)  and  antenna  2  (bottom) 


-O- . 


0  0.001  0.002  0.003  0.004  0.005  0.006  0.007  0.008  0.009  0.01 

Time  (s) 


<  ° 


0  0.001  0.002  0.003  0.004  0.005  0.006  0.007  0.008  0.009  0.01 

o 

g  Time  (s) 

o 

Fig.  37  Phase  calculated  from  filtered  IQ  data  (top),  phase  difference  (middle),  and 
corrected  phase  difference  (bottom) 
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Fig.  38  Four  AoA  trials  at  2  transmit  powers  (top)  and  average  error  vs.  transmit  power 
(bottom) 


In  the  test  setup,  the  transmit  antenna  gain  was  9  dB,  the  cable  loss  from  the  signal 
generator  to  the  antenna  was  5  dB,  and  the  distance  from  the  transmit  antenna  to 
the  receive  antenna  was  7.4  m.  Using  the  equation  for  free-space  path  loss, 

FSPL  =  20  log10  d  +  20  log10  /  +  20  log10  y,  (30) 

7.4  m  corresponds  to  a  path  loss  of  57.4  dB  at  2.4  GHz.  Including  the  other  gains 
and  losses  gives  a  total  loss  to  the  receive  antenna  of  -53.4  dB,  resulting  in  a  power 
at  the  receive  antenna  Prx  —  Ptx  —  53.4.  In  a  field  setting,  a  standard  dipole 
transmit  antenna  with  a  gain  of  2  dB,  no  cable  loss,  and  a  maximum  transmit  power 
of  10  dBm  would  give  a  total  transmit  power  of  12  dBm.  In  this  setting,  the  average 
error  versus  transmit  power  Ptx  in  Fig.  38  can  be  represented  as  distances  using 
Eq.  30  as 

d  —  10G2-ptx+53.4)/20+log10/+log10(47r/c)^  Q  ] ) 

shown  in  Fig.  39.  The  free-space  path  model  is  idealized,  and  might  lead  to  overly 
optimistic  range  calculations,  but  in  principle  Fig.  39  shows  that  AoA 
measurements  should  work  well  even  at  great  distances. 
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10 


Fig.  39  Average  AoA  error  vs.  distance 

The  third  experiment  studied  the  effect  of  antenna  spacing  on  AoA  error.  From 
Eq.  3,  we  see  that  as  the  antenna  spacing  increases,  the  AoA  becomes  less  sensitive 
to  variations  in  A <p.  Therefore,  greater  antenna  spacings  should  produce  less  AoA 
error.  This  theory  is  confirmed  by  the  bottom  plot  of  Fig.  40,  showing  AoA  error 
decreasing  as  the  antenna  spacing  increases.  The  top  plot  shows  4  trials,  2  at 
6.65-cm  spacing  and  2  at  38.6  cm.  While  the  smaller  spacing  produces  more  error, 
the  error  is  repeatable.  Thus,  even  constrained  applications  that  must  use  smaller 
antenna  spacings  may  be  able  to  achieve  accurate  measurements  if  the  error  can  be 
characterized  and  removed. 
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Fig.  40  Four  AoA  trials  using  2  antenna  spacings  (top)  and  average  error  vs.  antenna 
spacing  (bottom) 
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In  all  of  the  experiments  described,  the  unknown  SDR  phase  bias,  A(pbias,  was 
removed  in  postprocessing  by  subtracting  the  A <p  measurement  at  6  —  0°,  denoted 
as  A0O,  from  all  of  the  phase  measurements.  In  practice,  however,  this  bias 
information  is  unknown  and  requires  a  calibration  method  to  be  removed.  One 
method  is  to  transmit  a  signal  at  a  known  location  to  remove  the  bias.  A  third 
antenna  was  placed  in  the  center  of  the  antenna  array,  as  shown  in  Fig.  15.  A 
transmit  port  on  the  E310  was  used  to  transmit  from  this  third  antenna  and  receive 
on  the  other  2  antennas  using  a  calibration  program,  cal_test.py,  included  in 
Appendix  B.  This  A <fi  calibration  measurement  is  denoted  as  A<pcai.  A  normal 
measurement  was  then  performed  using  the  RF  signal  generator  source  at  9  —  0° 
to  find  A0O,  and  A<pbias  was  calculated  as 

Pbias  —  ^0  o  —  A<p  cai .  (32) 

At  a  later  date,  new  measurements  A <p'0  and  A(p'cal  were  taken.  Using  A(pbias  from 
Eq.  32,  the  new  A <p'0  measurement  was  estimated  as 

A0o  =  A(pbias  +  A  (p'cal.  (33) 

The  error  between  A(p'{)  and  A <p'0  was  0.04  radians,  corresponding  to  an  AoA  error 
of  0.17°  using  Eq.  2  and  the  standard  parameters  listed  in  Table  1.  This 
demonstrates  that  A(pbias  can  be  determined  and  removed  with  low  error. 

5.  Conclusion 


RDF  is  a  promising  technology  for  swarm  localization.  Initial  testing  in  an  anechoic 
chamber  with  an  array  of  2  antennas  showed  that  results  matched  well  with 
theoretical  calculations  over  a  range  of  antenna  spacings.  More  importantly, 
abnormities  were  repeatable,  showing  that  higher  accuracy  may  be  achieved 
through  calibration.  In  this  unique  application  where  the  swarm  has  control  over 
the  transmit  frequency,  multiple  frequencies  were  used  to  resolve  phase 
ambiguities.  The  use  of  a  narrowband  signals  allowed  filtering  to  remove  noise  and 
permitted  extremely  weak  signals  to  be  detected  and  processed.  A  calibration 
method  to  eliminate  the  unknown  phase  bias  was  demonstrated. 

Future  work  will  examine  the  use  of  additional  antennas  to  allow  for  higher  RDF 
resolution  and  the  determination  of  elevation  as  well  as  azimuth  angles.  High-speed 
RF  switches  can  be  used  to  sample  multiple  antennas  using  one  SDR  receiver 
port,22  allowing  a  single  E310  with  2  receivers  to  support  more  than  2  antennas. 
Subspace  methods  lend  themselves  to  analysis  of  the  multiple  signals  that 
additional  antennas  will  provide.  These  methods  can  be  explored  to  increase  RDF 
accuracy  and  determine  the  AoA  of  multiple  sources  simultaneously.23 
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Additionally,  experiments  in  realistic  environments  should  be  performed  and 
compared  with  results  from  the  anechoic  chamber.  Further  work  should  also 
consider  the  whole  system,  designing  a  scheme  for  the  relative  localization  of  an 
entire  swarm  based  on  RDF.  This  may  include  the  use  of  RF  sources  for  ranging 
and  communications  as  well  as  RDF.  Thus,  RDF  of  the  modulated  signals  used  in 
communications  should  also  be  researched.  System  design  will  need  to  consider 
medium  access  issues  and  effects  of  noise  and  jamming.  Localization  should  be 
modeled  together  with  navigation  to  characterize  swarm  performance  in  relation  to 
RDF  design  parameters.  Although  there  is  still  much  research  to  perform,  this 
report  has  taken  the  first  steps  toward  a  practical  RDF  system  for  swarm 
localization. 
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Appendix  A.  GNU  Radio  Radio-Direction-Finding  (RDF)  Program 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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#!/usr/bin/env  python2 
#  coding:  utf-8 

################################################## 

from  gnuradio  import  blocks 

from  gnuradio  import  eng_notation 

from  gnuradio  import  gr 

from  gnuradio  import  uhd 

from  gnuradio. eng_option  import  eng_option 

from  gnuradio. filter  import  firdes 

from  optparse  import  OptionParser 

import  time 

class  test_udp(gr.top_bk>ck): 
def _ init _ (self): 

gr.top_block. _ init  (self.  "test_udp") 

################################################## 

#  Blocks 

################################################## 
self.uhd_usrp_source_0  =  uhd.usrp_source( 

",".join((",  "")), 

cpu_format= "  fc  3  2 " , 
ot  w_format='  sc  1 6 
channels=range(2) , 

), 

) 

self.uhd_usrp_source_0.set_samp_rate(50e3) 
self.uhd_usrp_source_0.set_center_freq(2400e6,  0) 
self.uhd_usrp_source_0.set_gain(70,  0) 
self.uhd_usrp_source_0.set_center_freq(2400e6,  1) 
self.uhd_usrp_source_0.set_gain(70,  1) 

self  .blocks_udp_sink_0_0  =  blocks  .udp_sink(gr .  sizeof_gr_complex*  1 , 

'192.168.0.111',  15,  20000,  True) 

self.blocks_udp_sink_0  =  blocks.udp_sink(gr.sizeof_gr_complex*l, 
'192.168.0.111',  14,  20000,  True) 

################################################## 

#  Connections 

Approved  for  public  release;  distribution  is  unlimited. 

36 


################################################## 
self.connect((self.uhd_usrp_source_0,  0),  (self.blocks_udp_sink_0,  0)) 
self.connect((self.uhd_usrp_source_0,  1),  (self.blocks_udp_sink_0_0,  0)) 

def  main(top_block_cls=test_udp,  options=None): 

tb  =  top_block_cls() 

tb.start() 

try: 

raw_input('Press  Enter  to  quit:  ') 
except  EOFError: 

pass 

tb.stop() 

tb.wait() 

if _ name _ ==  ' _ main _ ': 

main() 
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Intentionally  left  blank. 
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Appendix  B.  GNU  Radio  Calibration  Program 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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#!/usr/bin/env  python2 
#  coding:  utf-8 

################################################## 

from  gnuradio  import  analog 

from  gnuradio  import  blocks 

from  gnuradio  import  eng_notation 

from  gnuradio  import  gr 

from  gnuradio  import  uhd 

from  gnuradio. eng_option  import  eng_option 

from  gnuradio. filter  import  firdes 

from  optparse  import  OptionParser 

import  time 

class  test_udp(gr.top_block): 
def _ init _ (self): 

gr.top_bk>ck. _ init _ (self,  "test_udp") 

################################################## 

#  Blocks 

################################################## 
self.uhd_usrp_source_0  =  uhd.usrp_source( 

",".join((", 

uhd.stream_args( 

cpu_format= "  fc  3  2 " , 
ot  w_format='  sc  1 6 
channels=range(2) , 

), 

) 

self.uhd_usrp_source_0.set_samp_rate(50e3) 
self.uhd_usrp_source_0.set_center_freq(2400e6,  0) 
self.uhd_usrp_source_0.set_gain(70,  0) 
self.uhd_usrp_source_0.set_center_freq(2400e6,  1) 
self.uhd_usrp_source_0.set_gain(70,  1) 

self.blocks_udp_sink_0_0  =  blocks.udp_sink(gr.sizeof_gr_complex*l, 
'192.168.0.111',  15,  20000,  True) 

self.blocks_udp_sink_0  =  blocks.udp_sink(gr.sizeof_gr_complex*l, 
'192.168.0.111',  14,  20000,  True) 

self.uhd_usrp_sink_0  =  uhd.usrp_sink( 
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uhd.stream_args( 
cpu_format="fc32", 
ot  w_format='  sc  1 6 
channels=range(  1 ) , 

), 

) 

self.uhd_usrp_sink_0.set_samp_rate(50e3) 

self.uhd_usrp_sink_0.set_center_freq(2400002000,0) 

self.uhd_usrp_sink_0.set_gain(70,0) 

self.analog_const_source_x_0  =  analog.  sig_source_c(0, 

analog . GR_C ONS T_W A VE,  0,  0,  0.7) 

################################################## 

#  Connections 

################################################## 
self.connect((self.uhd_usrp_source_0,  0),  (self.blocks_udp_sink_0,  0)) 
self.connect((self.uhd_usrp_source_0,  1),  (self.blocks_udp_sink_0_0,  0)) 
self.connect((self.analog_const_source_x_0,  0),  (self.uhd_usrp_sink_0,  0)) 
#self.connect((self.analog_const_source_x_l,  0),  (self.uhd_usrp_sink_0,  1)) 

def  main(top_block_cls=test_udp,  options=None): 

tb  =  top_block_cls() 

tb.start() 

try: 

raw_input('Press  Enter  to  quit:  ') 
except  EOFError: 

pass 

tb.stop() 

tb.wait() 

if _ name _ ==  ' _ main _ ': 

main() 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


AoA 

angle  of  arrival 

ARL 

US  Army  Research  Laboratory 

DDC 

digitally  down  converted 

GPS 

Global  Positioning  System 

IQ 

Inphase/Quadrature 

PC 

personal  computer 

PL 

programmable  logic 

PS 

processing  system 

RDF 

radio  direction  finding 

RF 

radio  frequency 

RFIC 

RF  integrated  circuit 

SDR 

software-defined  radios 

UAV 

unmanned  aerial  vehicle 

UDP 

user  datagram  protocol 

USRP 

Universal  Software  Radio  Peripheral 
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